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FIELD OF THE INVENTION 

The object of the present invention is a method for forming a neural network 
module, intended for real-time simulation of the flow mode, at any point of a pipe, of a 
multiphase fluid stream comprising at least a liquid phase and at least a gas phase, that 
5 is best suited to operating conditions and to a set of fixed physical quantities. 

BACKGROUND OF THE INVENTION 

Transporting hydrocarbons from production sites to treating plants constitutes an 
important link in the petroleum chain. It is a delicate link because of the complex 
+; interactions between the phases forming the transported effluents. The basic objective 

10 for operators is to reach an optimum productivity under the best safety conditions. They 
r therefore have to control as best they can the velocity and the temperature so as to avoid 

fy unnecessary pressure drops, unwanted deposits and unsteady-state flows. The method 

O that is generally used consists in modelling in the best possible way the transportation of 

complex multiphase streams so as to provide at all times an image of the flows in the 
15 various parts of the production chain, taking into account the precise constitution of the 
effluent, the flow rates, the pressures and the flow modes. 

There are currently various software modules for simulating the transport of 
complex multiphase streams, allowing to design suitable production equipments at an 
early stage. 

20 Patents US-5,550,761, FR-2,756,044 (US-6,028,992) and FR-2,756,045 (TJS- 

5,960,187) filed by the applicant notably describe modelling modules forming the 
TACITE model known to the man skilled in the art, allowing to simulate the transport 

j 
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of complex multiphase streams on steady or transient flow and capable of taking into 
account instability phenomena that occur because of the irregular geometry of the 
formation crossed by the pipe or of the topography thereof, referred to by specialists as 
« terrain slugging » or « severe slugging ». 

5 The simulation modules are as complex as the modelled phenomena. Precision and 

performances can only be obtained after a relatively long modelling time, which is not 
really compatible with real-time management. That is the reason why these modelling 
modules cannot be used as they are for real-time management of the production. It 
therefore appears necessary to use modelling methods offering a good compromise 

1 0 between calculating speed and accuracy of results. 

Patent application FR-00/09,889 filed by the applicant describes a method intended 
for real-time estimation of the flow mode, at any point of a pipe having a structure that 
can be defined by a certain number of structure parameters, of a multiphase fluid stream 
defined by several physical quantities and comprising liquid and gas phases. According 
15 to this method, the flow mode is modelled : 

- by forming a non-linear neural network with an input layer having as many inputs as 
there are structure parameters and physical quantities necessary for good estimation of 
the output, an output layer with as many outputs as there are quantities necessary for 
estimation of the flow mode, and at least one intermediate layer, 

20 - by forming a learning base with predetermined tables connecting various values 
obtained for the output data to the corresponding values of the input data, and 

- by determining, by iterations, weighting factors of the activation function allowing 
to properly connect the values in the input and output data tables. 
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Output data of the neural network is preferably analysed so as to sort, among the 
values of the output data of the neural network, only the pertinent data to be taken into 
account for iterative determination of the weighting coefficients of the activation 
function. 

5 SUMMARY OF THE INVENTION 

The method according to the invention allows to form a module (hydrodynamic or 
thermodynamic for example) intended for real-time simulation of the flow mode, at any 
point of a pipe, of a multiphase fluid stream comprising at least a liquid phase and at 
least a gas phase, so that it is best suited to fixed operating conditions concerning a 
10 certain number of determined structure and physical parameters relative to the pipe, and 
a set of determined physical quantities (hydrodynamic or thermodynamic quantities for 
example), with fixed variation ranges for the parameters and the physical quantities. 

It comprises using a modelling system based on non-linear neural networks having 
each inputs for structure parameters and physical quantities, outputs where quantities 
15 necessary for estimation of the flow mode are available, and at least one intermediate 
layer. The neural networks are determined iteratively so as to adjust to the values of a 
learning base with predetermined tables connecting various values obtained for the 
output data to the corresponding values of the input data. 

The method is characterized in that a learning base suited to the imposed operating 
20 conditions is formed and optimized neural networks best adjusted to the imposed 
operating conditions are generated. 
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In the case, for example, where the module is to be integrated in a general 
multiphase flow simulation model, both hydrodynamic and thermodynamic, the model 
is used to form the learning base so as to select the set of physical quantities that is best 
suited to the model, as well as the variation ranges fixed for said parameters and said 
5 physical quantities, and the optimized neural networks that best adjust to the learning 
base formed are generated. 

DETAILED DESCRIPTION 

General points concerning the formation of the neural networks system 

We consider a circulation of multiphase fluids in a pipe with at least a liquid phase 
and at least a gas phase, and we try to form a neural system allowing, from a certain 
number of geometrical and physical input data relative to the pipe and of physical data 
relative to the fluids, to give instantly, for each section of the fluid stream, an estimation 
of the flow mode or, from a certain number of data linked with the pressure and 
temperature conditions of the fluid, and with the composition of the fluid, to give the 
thermodynamic behaviour of the fluid. The whole of this input data constitutes a 
database. 

1) Input and output data 

The input data is for example : 

- geometrical and physical data of the pipes : diameter, roughness and angle of 
20 inclination of the pipe, and thermodynamic data : density of the gas, density of the 
liquid, viscosity of the gas, viscosity of the liquid, etc. ; the whole of this data 
constitutes an n-uplet, and 
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- hydrodynamic data characterizing the mixture : gas/liquid surface tension, volume 
fraction of gas, barycentric velocity of the mixture, etc., 

or : 

- thermodynamic data relative to the environment of the fluid ; pressure, temperature, 
5 the whole of this data constituting an n-uplet, as well as data characterizing the 

composition of the mixture : molar mass, critical properties, enthalpy coefficients, etc. 

The model produces for example at the output the hydrodynamic behaviour of the 
effluents, and notably the flow regime. It evaluates and delivers at two main outputs 
hydrodynamic data in the part of the pipe where the flow type is to be determined, the 
10 velocity difference dV between gas and liquid for example, or the stratified type flow 
fraction P (P e [0 ;1]) (type where the liquid phase flows in the lower part of the pipe). 
Other quantities qualifying the flow type can be calculated from these two outputs. 

In applications to thermodynamics, the model produces at the output the 
thermodynamic quantities of the effluents such as, for example, the number and the 
15 nature of the phases in presence or the molar fractions of the phases. 

2) Structure of the system of networks 

The function to be modelled by the hydrodynamic or thermodynamic model 
involves various subfunctions that can be solved by using, if necessary, several 
networks playing a weU-determined part within the model. For example, a network can 
20 be dedicated to the non-linear function regression, while another network is dedicated to 
the classification of the various flow regimes. Besides, the connections to be established 
between the networks or between the module and its outer environment can require 
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complementary data processing (data normalization, denormalization, etc.). The system 
is therefore referred to as based on neural networks. 

In order to connect the whole of the hydrodynamic or thermodynamic input and 
output data, a system of neural networks is formed, these networks being preferably of 

5 MLP type, well-known to the man skilled in the art, since it is particularly well-suited 
for physical phenomena modelling. Its structure allows to describe the dynamic as well 
as the static components of the phenomena, even by fixing, if necessary, some of its 
parameters at a reified value, therefore physically representative. Thus, knowing 
physical equations that govern the flows allows for example to enrich the network and 

10 to best adapt it to the physical phenomena modelled thereby. 

The neural network comprises (Fig.l) three layers for example : the input layer 
consisting of ten neurons corresponding to the ten data (mentioned above) of the 
complete physical model, an output layer consisting of two neurons corresponding to 
the two parameters dV and (3 sought, and an intermediate layer, referred to as hidden 
15 layer, whose number of neurons N c is optimized. The network is totally connected. The 
non-linearity of this network is obtained by a sigmoid activation function governing the 
behaviour of the neurons in the hidden layer. The neurons of the output layer can be 
selected linear or non-linear. The activation function can be the identity function for 
example. 

20 3) Learning 

The weights of this structure are determined at the end of a learning stage ; during 
this stage, the network is supplied with a set of data forming the learning base of the 
network, and the structure and the weights of the network are optimized by minimizing 
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the errors observed for all the samples of the base, between the output data resulting 
from network calculations and the data expected at the output, given by the base. The 
errors can be the absolute errors between the input and output quantities or the relative 
errors, according to the performance desired for the network. 

5 The generalization powers of the network are then tested from its capacity to 

properly calculate the two outputs for inputs that are unknown thereto. 

Implementation particularities 

The system based on neural networks modelling the hydrodynamic module or the 
thermodynamic module is defined by a set of parameters that the learning base allows to 
10 determine and to fix. These parameters are adaptable to the constraints imposed by the 
use to which this module is dedicated. 

If the module is intended for use in a general context, the learning base used has to 
be exhaustive, i.e. include all the available data concerning the various inputs/outputs ; 
the module can thus represent all of the cases that can be encountered in the field. 

15 If the module is intended for use in a more specialized framework or context, 

production of an oil field for example, adjustment is carried out on a more restricted 
learning base containing only the data that can be encountered in this field alone : type 
or structure of the pipes or pipelines, topographic variations, more restricted ranges of 
variation of certain hydrodynamic or thermodynamic quantities, etc. The parameters of 

20 the neural networks are in this case calculated so as to obtain a product specific to a 
given configuration, in order to represent in the best possible way the flows in operating 
ranges specific to this field. 



8 



This module is thus adaptable to the various objectives fixing the use thereof. 
Environment of the neural network module 

The module developed can be used alone, without associated modules. In this case, 
it is adjusted to databases formed according to the desired use, according to whether the 
5 context is general or specific. It is independent and allows to estimate the hydrodynamic 
or thermodynamic quantities likely to allow monitoring of the behaviour of effluent 
flows in pipelines, from a learning base containing physical and thermodynamic 
quantities data provided otherwise. 

O 

The module can also be implemented in any model (set of simulation software 
10 modules capable of simulating the transient and/or steady behaviour of multiphase 
flows, such as for example the TACITE model mentioned above). In this case, the 
hj module is adjusted to reference models representing the current state of the art, and to 

M= associated databases that widely cover the whole of the definition domain of the n-uplet 

formed by the input data. This adjustment thus takes into account the specificities of the 
15 software environment in which it is included. 

To form the learning base that is best suited to the model where the neural system is 
to be integrated, the model is advantageously used to connect the input and output data. 

Examples of use 

The module based on neural networks is suited for any application requiring 
20 hydrodynamic representation of effluents : formation of flow predictive simulation 
codes, real-time simulation or training simulation, etc. Active production control 
methods can also be based on such a module. 
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Optimization examples 

Consider for example the case where a network giving the classification of the flow 
type of multiphase fluids is to be optimized : stratified type or class CI flow (defined 
above), intermittent type or class C2 flow (characterized by a succession of liquid and 
5 gaseous slugs), or dispersed type or class C3 flow (the liquid being carried along in 
form of fine droplets), according to inputs as described in paragraph « General points 
concerning the formation of the network ». 

A neural network is first formed by means of a general base. The resulting optimum 
structure obtained contains 30 neurons, with a single hidden layer. The network 
10 confusion table obtained (Fig.2) expresses the classification percentage obtained in the 
three output classes of the network for each class of the database. 

The database used is then restricted to some particular configurations, which reduce 
the ranges of variation of the various inputs of the network. The optimum neural 
network consists in this case of 28 neurons and comprises a single hidden layer. The 
15 classification result is shown by the confusion table in Fig.3. 

It can immediately be seen that the good classification results distributed over the 
diagonal are greatly improved in relation to the previous case, formed by means of a 
less selective base, less suited to the configuration to which the module is applied. 

Adjustment of the database to the conditions of use thus allows, on the one hand, to 
20 modify the structure of the network and, on the other hand, to optimize the 
generalization results obtained. 



